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ABSTRACT 

We report the discovery of HATS-lb, a transiting extrasolar planet orbiting the moderately bright 
V=12.05 G dwarf star GSC 6652-00186, and the first planet discovered by HATSouth, a global network 
of autonomous wide-field telescopes. HATS-lb has a period of P ss 3.4465 d, mass of M p « 1.86 Mj, 
and radius of R. p w 1.30 The host star has a mass of 0.99 M©, and radius of 1.04 Rq. The discovery 
light curve of HATS-lb has near continuous coverage over several multi-day periods, demonstrating 
the power of using a global network of telescopes to discover transiting planets. 

Subject headings: planetary systems — stars: individual (HATS-1, GSC 6652-00186) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

The detection and study of extrasolar planets has be- 
come one of the fastest developing fields in astrophysics. 
This has been fueled by the extraordinarily rapid rate of 
exoplanetary discovery. The two most productive discov- 
ery techniques have been the radial velocity (RV) method 
(measuring the reflex RV signature of a planetary orbit 
in its parent star) and the transit method (measuring 
the slight decrease in stellar brightness that occurs if the 
planet happens to pass across its parent star's disk as 
viewed from the Earth). 

Even though among confirmed planets, much more 
were detected by the RV method than the transit 
method, the latter enables many studies that are not 
possible for RV detected planets, since for each transit- 
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ing system a lot more information is available than for an 
RV only system, like the planet radius, an unambiguous 
mass and the inclination of the orbit relative to the line 
of sight. 

In addition, observations and measurements are pos- 
sible for transiting planets that cannot be done with 
RV-only planets. It is possible to measure the misalign- 
ment between stellar spin and planetary o rbits through 
the Rossiter-McLaughl i n (RM) effect (e.g.. iQueloz et al.1 
2000 ) IWinn et al.H2005l. [20061 [20071 12008^120091 12010al l 
20lTt INarita et al.l 120071 I2008L l2009bl Ia L l2010bMal ~and 
many others). Information can be obtained about 
the planetary atmosphere by measuring the depen- 
dence of the transit depth ( plane t radius) on wave l ength 
/e.g.. iCharbonneau et all l200l iRedfield et all l200l 



ISnellen et alJl2008HSing et alJl2008n "2009. 2011a.b). The 
surface brightness and atmospheric temperature of the 
planet can be determined from observations of the sec- 
ondary eclipse, when the planet goes behind the star, 
and from the out of eclipse variat ions in the com - 
bined planet plus star brigh t ness Ce.g.. [Snellen fe Covin* 
20071 ISneTlen et all [20Ta IDeming et al.l I2007L "12011 



Charbonneau et alj |2008t iDemorv et al.l 12012ft . as well 



as a wide array of other theoretical and observational 
studies. 

For these reasons, increasing the sample of known tran- 
siting extrasolar planets, and extending the range of 
systems for which transits are detectable is very valu- 
able. The HATSou th network, described in detail in 
iBakos et. al.l (|2012f) (to be submitted), is designed to 
achieve both of these goals. With telescopes situated 
at three sites in the southern hemisphere, roughly 120° 
apart in longitude, HATSouth is capable of detecting 
transits at longer periods than any other ground based 
ne twork present l y in o peration. In addition, as detailed 
in IBakos et. al.l (|2012f ) , it is one of the most sensitive 
wide-field, ground based searches to small planets. Since 
a large part of the sky is being surveyed, many bright 
stars will be checked for transits. As a consequence, the 
HATSouth survey is expected to find many of its planets 
around bright stars, which makes the detailed follow up 
studies mentioned above feasible. 
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In this paper we present the first planet to come from 
the HATSouth survey: HATS-lb. Since this is the first 
transiting planet discovered by th esurvey, we describe 
in detail the data analysis methods and the procedures 
used to confirm the planetary nature of the object and 
derive the stellar and planetary parameters. 

The layout of the paper is as follows. In Section [5] we 
report the detection of the photometric signal and the 
follow-up spectroscopic and photometric observations of 
the host star, HATS-1. In Section [3] we describe the 
analysis of the data, beginning with the determination 
of the stellar parameters, continuing with a discussion 
of the methods used to rule out nonplanetary, false pos- 
itive scenarios which could mimic the photometric and 
spectroscopic observations, and finishing with a descrip- 
tion of our global modeling of the photometry and radial 
velocities. Our findings are discussed in Sectional 

2. OBSERVATIONS 
2.1. Photometric detection 
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Fig. 1. — Unbinned instrumental r-band light curve of HATS- 
1 folded with the period P = 3.4464588 days resulting from the 
global fit described in Sec tion [3] The solid line shows the best-fit 
transit model (see Section l3.3l l. In the lower panel we zoom-in on 
the transit; the dark filled points here show the light curve binned 
in phase using a bin-size of 0.002. 

Table [T] summarizes the HATSouth discovery observa- 
tions of HATS-lb. Observations were obtained over an 
eight month period from 2010 Jan through 2010 Aug 
using the HS1 (HATSouth 1) instrument in Chile, the 
HS3 instrument in Namibia, and the HS5 instrument 
in Australia. The observations, reduction, and analy- 
sis of_ttie_ljATSouyi_data were carried out as described 
in iBakos et. al.l (|2012f) . We detected a significant tran- 
sit signal in the light curve of GSC 6652-00186 (also 
known as 2MASS 11420608-2321174; a = ll h 42 m 06.12s, 
5 = -23° 21'17.4"; J2000; V= 12.053 ± 0.017mag APASS 
DR5, e.g. lHenden et alj|2009t see Figure p. 

2.2. Spectroscopy 

Table [2] summarizes the follow-up spectroscopic obser- 
vations which we obtained for HATS-1. Initial medium 
and low-resolution "reconnaissance" observations were 
obtained using th e Wide Field Spectrograph (WiFeS; 
iDopita et aTll2007D on the ANU 2.3 m telescope. These 
observations were used to rule out various common false 
positive scenarios (blends between an eclipsing binary 
and a giant star, or F-M type binary systems). We 



then obtained higher resolution, and higher RV preci- 
sion observations with the Coralie spectrograph on the 
Euler 1.2 m telescope at La Silla Observatory, FEROS 
on the MPG/ESO 2 2 m telescope at La Silla, FIES 
(|Diupvik fc And ersen 2010) on the 2.5 m Nordic Optical 
Telescope, and CYCLOPS on the 3.9 m Anglo- Australian 
Telescope (AAT) at Siding Spring Observatory (SSO). 
Below we describe the observations and reduction tech- 
niques for each spectroscopic instrument used. 

2.2.1. ANV 2.3m/ WiFeS 

High signal to noise (S/N), medium resolution recon- 
naissance spectroscopy was performed with WiFeS on 
the ANU 2.3 m telescope. WiFeS is an image slicing in- 
tegral field spectrograph with 25 slitlets each sampled 
at 0.5". Our observations, made with the red arm of 
the spectrograph, employed the R7000 grating and the 
RT480 dichroic, giving a spectral resolution of A/AA = 
7 000 and velocity dispersion of 21.6kms~ 1 pixel~ 1 in 
the wavelength region 5 200 - 7 000 A. The 4K x 4K 
CCD was read out in 2x binning in the spatial direction, 
with a detector gain of 0.9 e~ ADU -1 , and read noise of 
5e"s _1 pixel -1 . Only half the CCD containing the stel- 
lar signal was read out to reduce the readout time. RV 
standard star exposures were taken on each night during 
twilight. 

Reduction of the spectra was performed using the 
IRAF package CCDPROC, and extracted in TWOD- 
SPEC. Spectra from the three brightest image slices 
were extracted individually. Wavelength calibrations 
were performed with bracketing Ne-Ar arc lamp expo- 
sures. Atmospheric Oxygen B band lines in the region 
6 882 — 6 906 A were used to provide a first order wave- 
length correction. Object spectra were cross correlated 
against RV standard star spectra in the wavelength re- 
gion 5 270 — 5 650 A. Velocity measurements were ob- 
tained for each permutation of image slice and standard 
star spectra cross correlation, performed using the task 
FXCOR. The template and object spectra were Fourier 
filtered to remove low and high frequency noise, and a 
Gaussian fit was performed on the peak of the cross corre- 
lation function to derive the velocity shift. The final RV 
measurements were calculated as the mean of these veloc- 
ities weighted according to the cross correlation function 
peak and object spectrum S/N. 

The RV root mean square (RMS) scatter of multi- 
epoch WiFeS observations is 1 kms -1 , as measured from 
RV standard stars. We scheduled our observations at 
orbital phase quadrature of the candidate, where the ve- 
locity difference is the greatest. The results from five 
observations (with mean exposure time of 400 s) show no 
RV variation greater than 2kms _1 , indicating that the 
candidate was not an unblended eclipsing binary. 

We also obtained a single (700 s exposure) low res- 
olution (A/AA = 3 000) spectrum which was used to 
estimate T c ff* and logg*. We use the blue arm of 
WiFeS with the B3000 grating and RT560 dichroic, 
which gives a wavelength coverage of 3500 - 6000 A at 
A/AA = 3 000 resolution. Spe ctrophotometric standard 
stars from Ha muv et al.l (|1994f) are observed for flux cali- 
bration, and a Ne-Ar lamp is used for wavelength calibra- 
tion. Reductions are performed using the IRAF packages 
CCDPROC, KPNOSLIT and ONEDSPEC. The four 



HATS- lb 



3 



TABLE 1 

Summary of photometric observations 



Facility 


Date(s) 


Number of Images a Cadence (s) b 


Filter 


HS-1 


2010 Jan-2010 


Aug 


6018 


276 


Sloan 


r 


HS-3 


2010 Jan-2010 


Aug 


5391 


281 


Sloan 


r 


HS-5 


2010 Jan-2010 


Aug 


655 


271 


Sloan 


r 


Swope/SITe3 


2011 May 25 




57 


196 


Sloan 


i 


FTS/Spectral 


2011 May 28 




68 


60 


Sloan 


r 


FTS/Spectral 


2011 Jul 05 




71 


64 


Sloan 


i 


MPG/ES02.2/GROND 


2012 Jan 21 




78 


148 


Sloan 


9 


MPG/ES02.2/GROND 


2012 Jan 21 




80 


148 


Sloan 


r 


MPG/ESQ2.2/GROND 


2012 Jan 21 




78 


148 


Sloan 


i 



a Excludes images which were rejected as significant outliers in the fitting procedure. 

The mode time difference between consecutive points in each light curve. Due to visibility, weather, 
pauses for focusing, etc., none of the light curves have perfectly uniform time sampling. 



brightest image slices are summed, and the object spec- 
trum is divided by a black body white dwarf spectrum to 
remove the instrument sensitivity function. Flux calibra- 
ti on is pe r forme d according to the methodology set out 
in iBessell (|1999f ). The spectrum is ma tched to a grid of 
synth etic MARCS model atmospheres (jGustafsson et al.l 
I2008D . fitting for T e g+, log g*, [Fe/H], and interstellar red- 
dening (E(B - V)), at 250K, 0.5dex, and 0.01 magni- 
tude intervals respectively. Regions sensitive to changes 
in logg, such as the B aimer jump, the MgH feature at 
4 800 A, and the Mg b triplet at 5 170 A, are weighted 
preferentially. Based on the approximate log g+ and T c ff*, 
we find that HATS-1 is a dwarf. 

2.2.2. Euler l.Sm/Coralie 

We used the CORALIE spectrograph mounted on the 
Swiss 1.2 m Leonard Euler telescope in the ESO La Silla 
Observatory to measure 6 radial velocities (RVs) for 
HATS-lb during a run on 2011 May 19-20. CORALIE is 
a fiber- fed Echelle spectrograph with a design similar to 
ELODIE (Bar anne et al.l 1 19961) an d originally installed 
on April 1998 (jQueloz et al.l 120011) . It was refurbished 
in 2007 to increase its optical efficiency; the refurbished 
instrument provides a resolution of ~ 60 000. Spectra 
of our targets are taken in simultaneous Thorium- Argon 
(ThAr) mode, whereby one of the fibers observes the tar- 
get and the other a ThAr lamp (hereafter the object and 
comparison fibers, respectively). 

An automated pipeline was developed at Pontificia 
Universidad Catolica de Chile to extract and calibrate 
the CORALIE spectra, following in lar ge part the re- 
ductio n scheme for ELODIE described in Baranne et al. 
(1996). We provide a brief description of the reduction 
procedure here, full details will be given elsewhere. A 
night of observations include the following sets of im- 
ages: (i) science frames, where the object fiber observes 
the target and the comparison fiber is illuminated by the 
ThAr lamp; (ii) ThAr frames, where both fibers are il- 
luminated by the ThAr lamp; (iii) flats, where each of 
the fibers in turn is illuminated by a lamp at the begin- 
ning of the night. The pipeline trims and bias subtracts 
all frames. Flats are then median combined and used to 
find and trace each of the Echelle orders; the traces are 
fit by polynomials of order 4 and stored for extraction 
of the ThAr and science frames. We extract 70 orders 
from the object fiber and 51 from the comparison one, 
the wavelength coverage of the object fiber is «3 850- 



6 900 A. All ThAr frames are then extracted using the 
traces obtained from the flats, ThAr lines are identified 
and then fit by Gau ssians. Only lines from the list of 
ILovis fc Pepi (l2007f) are used. With the line centroids 
and wavelengths given in ILovis fc Pepd (|2007fl in hand, 
a global (i.e., for all orders simultaneously) wavelength 
solution of the form A = P(x,m)/m is then fit, for the 
science and comparison fibers independently, where x is 
pixel position of ThAr lines along the dispersion direc- 
tion, to is the Echelle order number and P(x, to) is a 
polynomial of order 3 in X and 5 in to. When fitting the 
wavelength solution we reject lines iteratively until the 
rms around the fit is acceptable. 

The spectrum in each Echelle order of the object fiber 
in the science frames i s extra c ted us ing the optimal ex- 
traction algorithm of iMarshl (|1989f ). The comparison 
ThAr spectra are extracted via a simple extraction and 
are used to calculate the velocity shift 6v of the compar- 
ison fiber of each science frame with respect to the com- 
parison fiber of the ThAr frame nearest in time. The 
wavelength calibration for the orders corresponding to 
the object fiber in the science frames is done by applying 
the wavelength solution of the object fiber of the ThAr 
frame nearest in time shifted by Sv. Finally, we apply 
a barycentric correction computed at the flux-weighted 
mean of the exposure and calculated using the JPL So- 
lar System ephemeris. With the spectra calibrated and 
corrected to the barycenter of the Solar System, a cross- 
correlation func tion is calculated usi ng a binary mask 
as described in iBaranne et al.l (|1996|) . A spectral typ- 
ing module within the pipeline gives estimates of T e ff, 
log(g), wsin(i) and [Fe/H] for each target, and we use 
these to choose an appropriate cross-correlation mask 
within a set of three choices: G2, K5, M2f TT l. The re- 
sulting CCF is fit with a Gaussian and the m ean re- 
turned by the fit is the measured RV. Following iQuelozl 
(1995) we estimate the uncertainty in the measured RVs 
via the scaling ctrv = a + [(1.6 + 0.2w)b/sn], where sn 
is the signal-to-noise ratio at 5 130 A , w is the width 
of the CCF measured in pixels, and a w 6ms~ 1 and 
b « 100 ms -1 are constants that are determined via sim- 
ulations and depend slightly on spectral type. We have 

11 The masks are the same as those used with the HARPS 
spectrograph, with the "transparent" regions of the masks made 
wider to match the lower resolution of CORALIE with respect to 
HARPS. Once a mask is chosen for a given target it is always 
maintained for further observations of it. 



4 



Penev et al. 



monitored the RV standard HD 72673 for « 2 years, ob- 
taining 40 RV measurements with S/N>60 per pixel at 
5 130 A. This star is known to have a long-term RV rms 
of ~ 2 ms" 1 from HIRES/Keck observations (Andrew 
Howard, private communication). Our CORALIE RV 
measurements have an RMS of 8 ms" 1 , which showcases 
the long-term velocity precision we can currently achieve 
with CORALIE for bright targets. The observations of 
HD 72673 show that even at high S /N we do not reach the 
level implied by the constant a in the expression for urv 
but rather we are limited to 0Rv,min = 8ms" 1 , which 
we adopt as our lower limit in the velocity uncertainty 
estimates. 

2.2.3. NO T 2.5 m/FIES 

We obtained two observations of HATS-1 using the 
Fibre-fed Echelle Spectrograph (FIES) on the 2.5 m 
Nordic Op tical Telescope (NOT) locat ed on the island of 
La Palma (jDiupvik fc Andersen|[20Tol) . The FIES obser- 
v ations were reduced to spectra following the procedure 
of lBuchhave et al.l (|2010l ). We used these spectra to mea- 
sure the atmospheric parameters of HATS-1, including 
the effective temperature T c g+, the metallicity [Fe/H], 
the projected rotational velocity wsini, and the stellar 
surface gravity logg*, by cross-correlating the observa- 
tions ag ainst a finely s ampled grid of synthetic spectra 
based on iKuruc z (2005) model atmospheres. The proce- 
dure, called Stellar Parameter Classification (SPC), will 
be described in detail in a fort hcoming paper (Bu chhave 
et al. , in preparation, see also IBakos et al.l 120121 for re- 
sults for other stars based on this method). 

2.2.4. ESO 2.2m/FEROS 

FEROS is the "Fiber-fed Extended Range Optical 
Spectrograph" located at the 2.2 m MPG/ESO telescope 
at the La Silla Observatory, Chile. The instrument covers 
the whole wavelength range between 3 500 and 9 200 A 
separated into 3 9 spectral orders at a re solving power 
A/ A A of 48 000 (jKaufer fc Pasauini|[l99l . making the 
analysis of the radial velocity and activity indicators such 
as Ca II H and K simultaneously possible (Sctiawan et al. 
l200l . 

FEROS is equipped with a double fiber system, pro- 
viding the opportunity to observe in object-sky or in 
object-calibration mode. The observations for HATS-1 
were taken in the object-calibration mode to guarantee 
the highest achievable precision with FEROS. In this 
mode, a ThAr+Ne lamp spectrum is recorded parallel to 
the target exposure, in order to calibrate the wavelength 
simultaneousl y. This method is mo deled on the ELODIE 
s pectrograph (|Baranne et al.1119961 ). 

Kau fer fc Pasquimi (|1998() report in the commissioning 
of FEROS at the 1.52 m telescope at La Silla observa- 
tory, Chile, a RV measurement e rror of ^23 ma" 1 . A n 
improved data reduction method (|Setiawan et alJ feOOO), 
reduced this error to ~10m/s, enabling FEROS to be 
used for planet searches. 

Between June 2011 and April 2012, a total of 21 spec- 
tra of HATS-1 were taken with FEROS during MPI 
guaranteed time (PI M. Mohler) and Chilean time (PI 
A. Jordan). 

The data reduction, including bias- and flatfield- 
correction, wavelength calibration and barycentric mo 
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Fig. 2. — Top panel: High-precision RV measurements for 
HATS-1 from CORALIE (dark filled circles), FEROS (open tri- 
angles) and CYCLOPS (filled triangles) shown as a function of 
orbital phase, together with our best-fit model. Zero phase cor- 
responds to the time of mid-transit. The center-of-mass velocity 
has been subtracted. Second panel: Velocity O — C residuals from 
the best fit. The error bars include a component from astrophysi- 
cal/instrumental jitter allowed to differ for the three instruments. 
Third panel: Bisector spans (BS), with the mean value subtracted. 
Note the different vertical scales of the panels. 

tion correction, was done using the MIDAS data reduc- 
tion pipeline at the telescope. 

The radial velocity of each spectrum was calculated 
by cross-correlating the object spectrum with a syn- 
thetic one, generated w ith the program SPE CTRUM 
(IGrav fc Corballvl Il99l using Kurucz models (jKurucd 
119931 ). In order to find the stellar parameters to synthe- 
size the fitting spectrum for cross-correlation, the object 
spectra were analyz ed using the program SM E (Spec- 
troscopy Made Easy, iValenti fc Piskunovl (|1996f) )■ Since 
the instrument efficiency of FEROS drops significantly 
below 4 000 A and above ^7 900 A, only the spectral 
range between these limits was used for cross-correlation. 
Due to the fact that Balmer lines, areas of telluric lines 
and strong emission features can falsify the result, these 
areas were carefully excluded in the analysis. In total, 26 
out of 39 orders were used. The RV was determined by 
fitting a Gaussian function to the final cross-correlation, 
with the center of the Gaussian fit representing the RV 
of the star. 

2.2.5. AAT 3.9m/CYCLOPS 

High-precision radial velocity observations were ob- 
tained with the CYCLOPS fiber-system feeding the 
UCLES spectrograph on the 3.9 m AAT at SSO, Aus- 
tralia. UCLES is a cross-dispersed Echelle spectrograph 
located at the coude focus with a 79 1 mm" 1 grating. CY- 
CLOPS is a Cassegrain fiber-based integral field unit 
used to feed UCLES. It consists of a 15 element fiber 
bundle which re-formats a ~ 3" diameter aperture into a 
pseudo-slit 0.63" wide and 15 elements long. It delivers 
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TABLE 2 

Summary of spectroscopic observations 



Telescope /Instrument 


Date Range 


Number of Observations 


Resolution 


ANU 2.3m/WiFeS 


2011 May 10—15 


5 


7 000 


ANU 2.3m/WiFeS 


2011 June 05 


1 


3 000 


AAT 3.9m/CYCLOPS 


2012 Jan 5—12 


6 


70 000 


Euler 1.2 m/Coralie 


2011 May 19—20 


6 


60 000 


NOT 2.5m/FIES 


2011 June 15—18 


2 


67 000 


MPG/ESO 2.2m/FEROS 


2011 June 08— 2012 April 18 


21 


48 000 



TABLE 3 

Relative radial velocities and bisector span measurements of 

HATS-1. 



BJD 


RV a 


o"RV b 


BS 


BS 


Phase 


Instrument 


(2 4^4 000+') 


[ 111 S j 


I'm s -1 ") 


(m s ) 










1701.46966 


64.78 


24.00 


30.4 


84.4 





.523 


Coralie 


1701.49218 


53.78 


25.00 


103.3 


87.4 





.530 


Coralie 


1701.67397 


132.78 


28.00 


3.7 


105.8 





.583 


Coralie 


1701.69858 


132.78 


29.00 


0.2 


140.7 





.590 


Coralie 


1702.46676 


229.78 


23.00 


-312.9 


87.2 





.813 


Coralie 


1702.49046 


209.78 


23.00 


2.7 


64.3 





.819 


Coralie 


1721.48200 


-190.09 


61.09 


-193.4 


257.4 





.330 


FEROS 


1722.56100 


246.23 


53.72 


97.1 


53.8 





.643 


FEROS 


1723.60400 


28.04 


59.75 


-10.1 


100.3 





.946 


FEROS 


1724.61100 


-409.37 


58.84 


-350.1 


98.5 





.238 


FEROS 


1726.63000 


157.47 


47.67 


-85.7 


139.2 





.824 


FEROS 


1728.53900 


-146.13 


56.44 


-263.5 


117.4 





.377 


FEROS 


1735.51300 


-121.06 


81.73 


-130.1 


514.8 





.401 


FEROS 


1736.57100 


286.48 


50.45 


30.3 


63.6 





.708 


FEROS 


1737.46400 


75.76 


50.47 


14.1 


70.5 





.967 


FEROS 


1738.46000 


-90.05 


57.78 


-46.1 


72.1 





.256 


FEROS 


1901.84000 


281.61 


34.46 


55.4 


129.7 





.661 


FEROS 


1911.84100 


79.18 


39.24 


236.8 


493.2 





.563 


FEROS 


1932.13480 


-2.80 


39.60 


-240.8 


57.8 





.451 


CYCLOPS 


1933.09769 d 


-104.80 


115.00 


-49.3 


279.6 





.731 


CYCLOPS 


1933.18812 


208.10 


33.62 


283.5 


59.9 





.757 


CYCLOPS 


1936.86700 


274.56 


48.03 


37.1 


49.0 





.824 


FEROS 


1938.16979 


-265.20 


21.28 


6.0 


66.2 





.202 


CYCLOPS 


1939.19224 


-22.00 


15.65 


1.6 


45.9 





.499 


CYCLOPS 


1943.22390 


271.60 


17.65 


-0.9 


22.0 





.669 


CYCLOPS 


1991.75100 


108.73 


56.55 


-19.4 


65.4 





.749 


FEROS 


2021.75000 


-33.57 


57.03 


105.1 


42.7 





.453 


FEROS 


2024.79300 


-50.14 


90.16 


76.1 


160.6 





.336 


FEROS 


2025.78200 


222.49 


46.75 


54.7 


64.5 





.623 


FEROS 


2029.73900 


163.37 


65.29 


-189.6 


192.5 





.772 


FEROS 


2031.54600 


-256.83 


46.88 


197.9 


161.3 





.296 


FEROS 


2032.70200 


81.66 


40.39 


-77.7 


130.8 





.631 


FEROS 


2035.65100 


27.54 


69.87 


35.9 


72.2 





.487 


FEROS 



a The zero-point of these velocities is arbitrary. An overall offset 7 re i fitted 
separately to the FIES, Keck and CYCLOPS velocities in Section l3.3l has been 
subtracted. 

^ Internal errors excluding the component of astrophysical/instrumental jitter 
considered in Scction |3.3l 



a spectral resolution of A/AA ^ 70 000 over 19 Echelle 
orders, with a total wavelength range of 4 540 - 7 340 A. 

The EEV2 2K x 4K CCD is readout with x2 bin- 
ning in the spatial direction and the normal readout 
speed, delivering 1.3 e~ ADU -1 gain and a read-noise of 
3.19e _ pix _1 . Exposure times vary between 1200s and 
2 400 s depending on the brightness of the target and the 
weather conditions. We take a ThAr exposure before 
and after each science exposure to monitor the stability 
of the spectrograph. 

Data reduction is performed using custom MATLAB 
routines, which trace each fiber and optimally extract 
each spectral order. Wavelength calibration is performed 



by linear interpolation between the ThAr exposures on 
either side of the exposure. Using the IRAF task FXCOR 
we cross-correlate each of the individual fiber spectra 
against corresponding fiber spectra of a very high signal- 
to-noise radial velocity standard star exposure. Final ve- 
locities are determined as the weighted mean of all fibers 
over all orders with 3 sigma clipping. Typically 12 fibers 
over 17 orders are used, with the two reddest orders too 
heavily contaminated with telluric lines to deliver pre- 
cise RV measurements. RV uncertainties are calculated 
as the standard deviation of the scatter of the fibers over 
all the orders used. 
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2.3. Photometric follow-up observations 



2011 May 25 (SWOPE) 




-0.15 -0.1 -0.05 0.05 0.1 0.15 

Time from transit center (days) 

Fig. 3. — Unbinncd instrumental Sloan g-, r-, and i-band transit 
light curves of HATS-1. The dates and instruments used for each 
event are indicated. The light curves have been detrended using 
the EPD and TFA processes. Curves after the first are shifted for 
clarity. Our best fit is shown by the solid lines. Residuals from 
the fits are displayed at the bottom, in the same order as the top 
curves. 

Follow-up photometric observations of HATS-1 were 
obtained using the SITe3 camera on the Swope 1.0 m 
telescope at LCO, the Spectral instrument on the 2.0 m 
Faulkes Telescope South (FTS) at SSO, and GROND on 
the MPG/ESO 2.2 m telescope at La Silla Observatory 
in Chile. Table [T] summarizes these observations. Below 
we describe the reduction procedures for each instrument 
used. 

2.3.1. Swope lm/SITe3 

The SITe#3 camera on the Swope 1 m telescope is a 
2048x3150 detector with a pixel scale of 0.435" and a 
field of view of 14.8' x 22.8'. Observations of a field cen- 
tered on HATS-1 were calibrated using standard methods 
and reduced to light curves follow ing a similar procedure 
to that used bv lBakos et al.l (|2010T l in the analysis of pho- 
tometric follow-up data for the HATNet pr oject . Briefly, 
we use a variation of the fitsh package ()Pall I2012T ) to 
identify sources on the images, cross-match them be- 
tween images, perform aperture photometry on the im- 
ages using a range of apertures, collect light curves for 
the detected sources, and perform an ensemble magni- 
tude fit to correct the light curves for variations that are 
common to multiple sources. Auxiliary parameters de- 
scribing the shape of the stellar profile, the position of 
each source on each image, the sky background, and the 
local variation in the sky background, are also recorded. 



These are used in the light curve fitting procedure (§ I3.3[) 
to correct for systematic variations in the target star that 
are not removed by the ensemble magnitude fit. 

2.3.2. FTS 2 m/ Spectral 

Two partial transits of HATS-lb were observed with 
the FTS, which is part of the Las Cumbres Global Tele- 
scope (LCOGT) Network. We used the 4Kx4K "Spec- 
tral" imaging camera, featuring 0.15" pixels and a 10' 
FoV. Exposures of 40 s were taken in both r-band (2011 
May 28) and i-band (2011 July 5), binning 2 x 2 to 
reduce readout time, and slightly de-focusing the tele- 
scope to avoid saturation. Raw data is reduced via the 
automated reduction pipeline provided by LCOGT for 
FTS data. Photometry is performed by an automated 
pipel ine based on Source Extractor (|Bertin k. Arnoutsl 
I1996T ) aperture photometry and calibrated using selected 
neighboring reference stars. Details of these observations 
are set out in Table 2J and photometry of the two transit 
ingresses observed is presented in Figure [T] 

2.3.3. ESO 2.2m/ 'GROND 

A complete transit of HATS-lb was observed on UT 
2012 January 21 in the Sloan g (A = 4 550 A), r (A = 
6 270 A) and i-band (A = 7 630 A) of the GROND instru- 
ment (iGreiner et al.ll2008f l. mounted on the MPG/ESO 
2.2m telescope at the ESO La Silla Observatory (Chile). 
The simultaneous multiband time scries was obtained 
applying a moderate telescope defocus to minimize flat- 
ficlding errors and to avoid saturation due to abrupt see- 
ing decrease. The field-of-view (FOV) of GROND in the 
three optical channels is 5.4' x 5.4' with a pixel scale of 
0.158" pixel -1 , which allowed monitoring of the target 
and a few nearby reference stars. We acquired data for 
nearly 3.5 hr time interval, covering the predicted transit 
of HATS-lb. 

The data reduction and analysis of the photometric 
data was performed using a customized pipeline, follow- 
ing standard procedures. The masterbias and master- 
flat frames were used to perform debias and flat-fielding 
of the science images. DAOPHOT aperture photometry 
was performed as implemented in the IDL APER func- 
tion to obtain light curves of HATS-1 and the reference 
stars in GROND's FOV with three aperture radii (27.5, 
30 and 32.5 pixels). Finally a differential photometry 
was performed relative to the reference stars. 

3. ANALYSIS 
3.1. Properties of the parent star 

As mentioned in Section [2.2.3l we determine the atmo- 
spheric parameters of the host star HATS-1, including 
the effective temperature T c g+, the metallicity [Fe/H], 
the projected rotational velocity usinz, and the stellar 
surface gravity logg*, by applying SPC to the FIES ob- 
servations of this star. This analysis yielded the following 
initial values and uncertainties: T B s* — 5870 ± 100 K, 
[Fe/H] = -0.06 ± 0.12 dex, wsini = 2.52 ± 0.5kms~\ 
and log 5* = 4.52 ± 0.15 (cgs). 

The same parameters were independently determined 
by applying SME to the FEROS observations. The 
values obtained were T off * = 5790 ± 120 K, [Fe/H] = 
—0.11 ± O.ldex, vsini — 4.8 ± 0.5, kms -1 , and logg* = 
4.53 ± 0.17 (cgs). These measurements, apart from vsini 
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TABLE 4 

Differential photometry of HATS-1 



BJD 
(2 400 000+) 


Mag a 


o-Mag Mag(orig) b Filter 


Instrument 


55287.81456 


-0.01671 


0.00299 r 


HS 


55308.49356 


0.00495 


0.00267 ■ • ■ r 


HS 


55325.72598 


-0.00006 


0.00260 ■ • ■ r 


HS 


55384.31621 


-0.00972 


0.00255 ■ ■ ■ r 


HS 


55294.70845 


0.00056 


0.00244 ■ • ■ r 


HS 


55263.69047 


0.00764 


0.00313 ■•■ r 


HS 


55239.56544 


0.00769 


0.00264 ■ ■ ■ r 


HS 


55301.60184 


0.00695 


0.00246 ■ • ■ r 


HS 


55277.47684 


0.00689 


0.00280 ■ • ■ r 


HS 


55287.81767 


0.00245 


0.00299 ■ • ■ r 


HS 


Note. — This table is available in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its form and content. 
a The out-of-transit level has been subtracted. For the HATSouth light curve (rows 



with "HS" in the Instrument column), these magnitudes have been detrended using 
the EPD and TFA procedures prior to fitting a transit model to the light curve. 
Primarily as a result of this detrending, but also due to blending from neighbors, the 
apparent HATSouth transit depth is ~ 79% that of the true depth in the Sloan r 
filter. For the follow-up light curves (rows with an Instrument other than "HS" ) 
these magnitudes have been detrended with the EPD and TFA procedures, carried 
out simultaneously with the transit fit (the transit shape is preserved in this process), 
k Raw magnitude values without application of the EPD and TFA procedures. This 
is only reported for the follow-up light curves. 

are consistent within the quoted error bars. For the re- 
mainder of the paper the SPC results are adopted. The 
apparent inconsistency in the projected stellar rotation 
velocity is likely due to the mis-estimation of the quoted 
errors. For such slowly rotating stars, the uncertain- 
ties in both measurements are significantly larger than 
0.5kms _1 . 

Following Sozzctti ct al. (2007) we use the stellar mean 
density p*, together with T e ft* and [Fe/H] from SPC to 
constrain the fundamental properties of the star (mass, 
radiu s, age and lu minosity) based on the Yonsei-Yale 
(Y2. lYi et al.ll200~iT ) stellar evolution models. The stellar 
mean density is closely related to the normalized semi- 
major axis a/i?+, and is determined by modeling the light 
curve 13. 3p . The light curve model also depends on the 
quadratic limb darkening c oefficients, wh ich we initially 
fix to the values from the iClaretl ()2004D law using the 
initial Teg*, [Fe/H], and logy* values. The comparison 
with the stellar evolution models provides a more precise 
estimate of log g+ than can be determined from the spec- 
trum. We we then fix log g+ to this more precise value in 
a second iteration of SPC to determine revised estimates 
of the other atmospheric parameters (T e ff*, [F e /H], and 
wsini). These new atmospheric parameters are used to 
determine new estimates of the limb darkening coeffi- 
cients, which are used in a second iteration of the light 
curve modeling. We find that the value of logg* con- 
verges after two iterations. 

We collect the final values for the observed and derived 
stellar parameters for HATS-1 in Table [SJ The inferred 
location of the star in a diagram of a/i?* versus Teg*, 
analogous to the classical H-R diagram, is shown in Fig- 
ure |U The stellar properties and their la and 2a con- 
fidence e l lipsoi ds are displayed against the backdrop of 
lYi et al.l {2001) isochrones for the measured metallicity 
of [Fc/H] = —0.1, and a range of ages. 

3.2. Excluding blend scenarios 



4.0 

5.0 
6.0 



7.0 




6500 6000 5500 

Effective temperature [K] 



Fig. 4. — Model isochrones from Yi et al. (2001) for the measured 
metallicity of HATS-1, [Fe/H]= -0.1, and ages of 0.2 Gyr, and 1 to 
13Gyr in 1 Gyr increments (left to right). The adopted values of 
T e ff* and a/i?* are shown together with their 1<t and 2cr confidence 
ellipsoids. The open triangle shows the values from the initial SPC 
iteration. 

Most of the obvious astrophysical false positive scenar- 
ios which might produce the detected transit of HATS-1 
are ruled out by the spectroscopic observations discussed 
in Section 12.21 Some scenarios involving diluted eclips- 
ing binary systems, or transiting planet systems diluted 
by an additional star, may still be consistent with these 
observations. To rule these scenarios out we conduct 
detailed modeling of the light curves foll owing the pro- 
cedure described in lHartman et all (|2011| ). 

From the light curves alone we are able to rule out hi- 
erarchical triple star systems with greater than 4a confi- 
dence, and blends between a foreground star and a back- 
ground eclipsing binary with > 2.5a confidence. More- 
over, the only nonplanetary blend scenarios which could 
plausibly fit the light curves (scenarios for which the re- 
jection confidence is less than 5cx) are scenarios in which 
the two brightest stars in the blended system are nearly 
equal in brightness; such scenarios would have been easily 
detected from the spectroscopic observations (the spec- 
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TABLE 5 
Stellar, parameters for HATS-1 



Parameter 


Value 


Source 


Spectroscopic properties 




v;.,r. (K) 


5870 ± 100 


SPC a 


[Fe/H] 


-0.06 ±0.12 


SPC 


v sin i (km s 1 ) . . . 


2.52 ± 0.5 


SPC 


Photometric properties 






v ijriagj 




APASS 


a ^magj 


i o 74fl + 066 


APASS 


y v liltx 6/ 


12.356 ± 0.052 


APASS 


r (mag) 


11.866 ±0.045 


APASS 


i (mag) 


11.748 ±0.042 


APASS 


J (mag) 


10.924 ±0.024 


2MASS 


H (mag) 


10.634 ± 0.022 


2MASS 


K s (mag) 


10.580 ±0.023 


2MASS 


Derived properties 






M* (M ) 


0.986 ± 0.054 


YY+a/_R*+SPC b 


/>'• I 


1.038±°;£g 


YY+a/R*+SPC 


log g± (cgs) 


4.40 ±0.08 


YY+a/ii*+SPC 


L* (L ) 


1 10+ ' 33 

1 - lu -0.18 


YY+a/R*+SPC 


M v (mag) 


4.73 ± 0.24 


YY+a/ii*+SPC 


M K (mag,ESO) . . 


3.21 ±0.22 


YY+a/ii*+SPC 


Age (Gyr) 


6.0 ±2.8 


YY+a/R*+SPC 


Distance (pc) 


303±H 


YY+a/ii*+SPC 



a SPC — "Spectroscopic Parameter Classification" procedure ap- 
plied to the high-resolution FIES spectra (Buchhave et al., in prepa- 
ration). These parameters rely primarily on SPC, but have a 
small dependence also on the iterative analysis incorporating the 
isochrone search and global modeling of the data, as described in 
the text. 

b YY+a//i,+SPC = Based on the YY isochrones lYi et al .1120011) . 
a/R* as a luminosity indicator, and the SPC results. 

trum would have either been obviously double-lined, or 
there would have been several kms -1 changes in the RV, 
BS, and/or FWHM; the latter is computed, and checked 
for variations, but not shown). 

While our analysis rules out nonplanetary explanations 
for the observations, we are unable to disprove the pos- 
sibility that the object is an unresolved stellar binary 
system (or chance alignment of two stars) with one com- 
ponent hosting a transiting planet. If this is the case the 
true mass and radius of the planet would be larger than 
what we infer. However, lacking any positive evidence for 
multiple stellar components in the system, for the anal- 
ysis conducted in Section 15751 we assume that HATS-1 
consists of a single star with a transiting planet. 

3.3. Global modeling of the data 

We conducted a joint modeling of the HATSouth pho- 
tometry, follow-up photom etry, and RV data following 
the procedure described by [Bakos et al.l (|2010D . Briefly, 
we mod eled the light curves using analytic formulae 
based on iMandel fc Agoll (|2002l ) to describe the physical 
variation together with EPD and TFA to describe instru- 
mental variations, and we modeled the RV data w ith a 
Keplerian orbit with the formalism of iPall ( 20091) . We 
used the Downhill Simplex Algorithm (e.g. IPress et al.l 



1992) to optimize the parameters, followed b y a Markov - 
Chain Monte Carlo analysis (MCMC, e.g. iForl 12001 
to determine their uncertainties and correlations. The 
resulting geometric parameters pertaining to the light 
curves and velocity curves, as well as the derived physi- 
cal planetary parameters, are listed in Table [6j 



Included in this table is the RV "jitter" which was 
added in quadrature to the formal RV errors for the 
FEROS observations, such that x 2 /dof = 1 for this in- 
strument. The origin of this noise is unknown, it may 
either be to due astrophysical noise intrinsic to the star 
(potentially including additional planets in the system), 
or instrumental errors not included in the formal uncer- 
tainties. 

We find that the planet has a mass of M p = 
1.855±^Mj and radius of Rp = 1.302±£Jg Rj, which 
are typical values for a hot Jupiter. We allow the eccen- 
tricity to vary in the fit so as to allow our uncertainty 
on this parameter to be included in the uncertainties of 
the other physical parameters of the system, however we 
find e = 0.120 ± 0.092, which is consistent with zero ec- 
centricity. Our 99% upper-limit on the eccentricity is 
e < 0.4. 

4. DISCUSSION 

The HATSouth global network of telescopes relies on 
combining observations from three stations located at 
three different sites in the southern hemisphere, with a 
close to optimal longitude separations. This makes it 
possible to continuously (or nearly so if the night is short) 
observe a given portion of the sky to search for transits 
of planets in front of their parent stars. 

In this paper we have presented the first planet discov- 
ered by the HATSouth network. This discovery demon- 
strates that we are indeed able to successfully combine 
observations from multiple telescopes and sites to detect 
planetary transit signals. 

In many respects this is a very typical transiting planet, 
with an orbital period of P « 3.4465 days (close to the 
mode for transiting planets detected from the ground), 
mass M p w 1.86 Mj (slightly higher than typical), radius 
R p « 1.30 Rj (close to the median for its mass), around 
a star that is an almost exact analogue of the Sun. 

Figure O demonstrates the power of using a network of 
telescopes to search for transit signals among stars. It 
shows a ten day period (starting 2010 February 15) from 
the HATSouth discovery light curve containing 1351 ob- 
servations of HATS-1, i.e. a small fraction of the total of 
more than 12 000 observations. Over this period, a signif- 
icant fraction of the observations come from each of the 
three HATSouth sites, and a transit is detected by each 
station at a comparable significance. Despite the short 
nights in the southern hemisphere at this time, near con- 
tinuous coverage is achieved and every single transit in 
this time window is observed by the HATSouth network. 
In contrast, a search operating from a single location 
would only be a ble to detect one of t hese transits. 

As discussed in I Bakos et. al.l (|2012t ) this increased duty 
cycle should result in many more planet detections, es- 
pecially at longer periods and for smaller planetary radii, 
whic h are arguabl y the m ost valuable. 

In iBakos et. al.l (|2012[ ) detailed simulations were car- 
ried out showing that the distribution of planets ex- 
pected to be produced by the HATSouth survey peaks 
at planet radii between 1 and 2 Jupiter radii and or- 
bits with periods between 3 and 5 days. Not surpris- 
ingly, the first planetary system discovered by the net- 
work falls within these ranges. However, as discussed in 
that paper, these predicted distributions predict signifi- 
cantly enhanced sensitivity to planets with much longer 
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TABLE 6 
Orbital and planetary parameters 



orbital periods and smaller radii compared to a single 
site survey. 



Parameter Value 

Light curve parameters 

P (days) 3.446459 ± 0.000004 

T c (BJD) a 2455241.28722 ± 0.00071 

T w (days) a 0.1007 ± 0.0012 

T12 = T 34 (days) a 0.0181 ± 0.0011 

a/R* 9.20 ±0.82 

C/R* h 23.98 ±0.16 

R p /R* 0.1288 ± 0.0020 

b^acosi/R-t, 0.633±oo23 

i(deg) 85.6+?;l 

Limb-darkening coefficients c 

a g (linear term) 0.5230 

b g (quadratic term) 0.2574 

a r 0.3387 

b r 0.3343 

a,i 0.2560 

bi 0.3331 

RV parameters 

K (ms^ 1 ) 254.1 ±33.1 

ecosu -0.064tg gg4 

esinw 0.097^0073 

e 0.120 ±0.092 

RV jitter (ms" 1 ) 11 61.5 

Planetary parameters 

M P (Mj) l-855t°;^ 

R p t-302+o Jgg 

C(M P , R p ) e 0.53 

p p (gem- 3 ) 1.03 ±0.25 

logg p (cgs) 3.43 ±0.07 

a(AU) 0.0444 ± 0.0008 

T cq (K) 1359±^ 

e f 0.127 ±0.014 

(F) (lO^rgs-^m" 2 ) s 7.70l?;g 

a T c : Reference epoch of mid transit that minimizes the correlation 
with the orbital period. BJD is calculated from UTC. T14: total 
transit duration, time between first to last contact; T3.2 — T34: 
ingress/egress time, time between first and second, or third and 
fourth contact. 

Reciprocal of the half duration of the transit used as a 
j ump parameter in our MCMC analysis in place of a/ .R* . It 
is related to a/ i?* by the expression £ / R* — a/ i?* (2tt(1 + 
esinu))/(P\/l - b' 2 y/l - e 2 ) iBakos et alj|20lch . 
c Values for a quadratic law given separately for the Sloan g, r, and 
i filters. These values were adopted from the tabulations by | Cl aret 
(2004) according to the spectroscopic (SPC) parameters listed in 
Tabled] 

" This jitter was added to the FEROS measurements only. Formal 
uncertainties for the Coralic measurements were not determined 
from the pipeline, instead the uncertainties were fixed such that 
X 2 /dof — 1 for these observations, so that, by definition, no jitter 
was required. 

e Correlation coefficient between the planetary mass M p and radius 
Ftp . 

The Safronov number is given by G - | (VWc/^orb) 2 — 
(a/Rp){M p /M*) (see[H arisen & Barman 2007). 
^ Incoming fiux per unit surface area, averaged over the orbit. 



Development of the HATSouth project was funded 
by NSF MRI grant NSF/AST-0723074, operations are 
supported by NASA grant NNX09AB29G, and follow- 
up observations receive partial support from grant 
NSF/AST-1108686. Followup observations with the 
ESO 2.2m/FEROS instrument were performed un- 
der MPI guaranteed time (P087.A-9014(A), P088.A- 
9008(A), P089.A-9008(A)) and Chilean time (P087.C- 
0508(A)). A.J. acknowledges support from Fondecyt 
project 1095213, Ministry of Economy ICM Nuclei P07- 
021-F and P10-022-F, Anillo ACT-086 and BASAL 
CATA PFB-06. V.S. acknowledges support form BASAL 
CATA PFB-06. M.R. acknowledges support from a 
FONDECYT postdoctoral fellowship. R.B. and N.E. ac- 
knowledge support from Fondecyt project 1095213. This 
work is based on observations made with ESO Tele- 
scopes at the La Silla Observatory under programme 
IDs P087.A-9014(A), P088.A-9008(A), P089.A-9008(A), 
P087.C-0508(A), 089.A-9006(A), and observations made 
with the Nordic Optical Telescope, operated on the is- 
land of La Palma jointly by Denmark, Finland, Ice- 
land, Norway, and Sweden, in the Spanish Observato- 
rio del Roque de los Muchachos of the Instituto de As- 
trofisica de Canarias. This paper also uses observa- 
tions obtained with facilities of the Las Cumbres Ob- 
servatory Global Telescope. Work at the Australian Na- 
tional University is supported by ARC Laureate Fellow- 
ship Grant FL0992131. We acknowledge the use of the 
AAVSO Photometric All-Sky Survey (APASS), funded 
by the Robert Martin Ayers Sciences Fund, and the SIM- 
BAD database, operated at CDS, Strasbourg, France. 
M. Mohler wants to thank N. Piskunov and his group in 
Uppsala for introducing her to Spectroscopy Made Easy 
(SME). Operations at the MPG/ESO 2.2m Telescope 
are jointly performed by the Max Planck Gesellschaft 
and the European Southern Observatory. The imaging 
system GROND has been built by the high-energy group 
of MPE in collaboration with the LSW Tautenburg and 
ESO. We thank Timo Anguita and Regis Lachaume for 
their technical assistance during the observations at the 
MPG/ESO 2.2 m Telescope. We thank Miguel Roth, 
Francesco Di Mille and Rodolfo Angeloni for allowing 
us to use the Swope telescope on May 24/25 2011. We 
are grateful to P.Sackett for her help in the early phase 
of the HATSouth project. 



REFERENCES 



Torres, G., Pal, A., et al. 2010, ApJ, 710, 1724 
Hartman, J. D., Torres, G., et al. 2012, arXiv: 



Bakos, G. A. 
Bakos, G. A. 

1201.0659 
Bakos et. al. 2012 

Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Bessell, M. S. 1999, PASP, 111, 1426 



Buchhave, L. A., Bakos, G. A., Hartman, J. D., et al. 2010, ApJ, 
720, 1118 

Charbonneau, D., Brown, T. M., Noyes, R. W., & Gilliland, R. L. 

2002, ApJ, 568, 377 
Charbonneau, D., Knutson, H. A., Barman, T., et al. 2008, ApJ, 

686, 1341 
Claret, A. 2004, A&A, 428, 1001 



10 



Penev et al. 




HJD-2455243 Phase 
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